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RESUME 

Nous  avons  developpe  des  la  fin  des  annees  70 
dans  notre  logiciel  ELFINI  des  techniques 
mathematiques  originales  d'identification  de 
modeles,  en  particulier  pour  les  modeles 
dynamiques  Elements  Finis  et  pour  le  recalage  des 
modeles  de  charges  aeroelastiques  stationnaires  sur 
les  essais  en  vol. 

Nous  nous  sommes  maintenant  attaques  au 
recalage  des  charges  aerodynamiques 
instationnaires  sur  les  essais  de  vibration  en  vol. 

Sur  le  plan  mathematique  nous  utilisons  toujours  la 
meme  approche,  qui  diflere  notablement  des 
methodes  classiques  de  moindres  carres 
(minimisation  d'une  "distance"  calcul-mesure  en 
fonction  des  parametres  d'ajustement).  Nous 
preferons  utiliser  une  approche  type  optimisation 
quadratique  en  minimisant  une  "distance"  entre  les 
parametres  d’ajustement  et  leurs  valeurs  theoriques 
ou  presumees,  en  formant  la  solution  a satisfaire  les 
inequations  de  reconstitution  des  mesures  par  le 
modele  k une  precision  donnee. 

Parmi  tous  les  avantages  de  cette  technique  le 
principal  est  de  s'affranchir  du  probleme  des 
parametres  mal  observables. 

Nous  montrons  2 applications  en  aeroelasticite  : 

- Le  recalage  des  champs  de  pression 

aerodynamiques  stationnaires  par  les  mesures 
en  vol  de  reponses  de  jauges  de  contraintes  en 
manoeuvre,  illustre  par  un  exemple  issu  de  la 
mise  au  point  des  charges  du  RAFALE. 

- Le  recalage  des  forces  aerodynamiques 

instationnaires  sur  la  mesure  en  soufflerie,  mais 
transposable  en  vol,  des  frequences  et 
amortissements  d'une  maquette  aeroelastique 
d'empennage  type  Airbus.  Nous  montrons  que 
le  calcul  de  la  vitesse  critique  de  Flutter  est 
recale  a partir  d'essais  a vitesse  bien  inferieure. 

Nous  concluons  en  presentant  les  perspectives  de 
developpement  de  la  methode. 


ABSTRACT 

Since  the  end  of  the  70ies  we  have  developed, 
within  the  frame  of  our  ELFINI  software,  original 
techniques  for  mathematical  model  identification, 
in  particular  for  calibration  of  dynamic  Finite 
Element  models  from  ground  vibration  tests  and  of 
steady  aeroelastic  loads  from  flight  tests. 


Now  we  have  tackled  calibration  of  unsteady 
airloads  from  flight  vibration  tests. 


Mathematically  speaking,  we  keep  the  same 
approach,  which  differs  notably  from  classical  least 
square  methods  (minimization  of  a calculation  - 
measurement  "distance"  in  function  of  calibration 
parameters).  We  prefer  to  use  a quadratic 
optimization  type  approach  with  the  minimization 
of  a "distance"  between  calibration  parameters  and 
their  theoretical  or  presumed  values,  constraining 
the  solution  to  satisfy  measurement  reconstitution 
by  the  model  at  a given  accuracy. 

Among  advantages  of  this  technique,  the  principal 
is  to  get  rid  of  ill-observable  parameters. 


We  describe  two  applications  : 

- Calibration  of  steady  aerodynamic  pressure 
fields  from  flight  measurements  of  strain-gage 
responses  in  maneuver,  illustrated  by  an  example 
coming  from  the  calibration  of  Rafale  airloads, 


- Calibration  of  unsteady  airloads  from  wind 
tunnel  measurements,  transposable  to  flight 
measurements,  of  frequencies  and  dampings  of  an 
aeroelastic  dynamic  model  of  an  Airbus  type 
stabilizer.  We  show  that  critical  flutter  speed  is 
identified  from  tests  at  much  lower  speed. 


As  a conclusion  we  present  the  future  prospects  of 
the  method. 


Paper  presented  at  the  RTO  AVT  Specialists’  Meeting  on  “Structural  Aspects  of  Flexible  Aircraft  Control”, 
held  in  Ottawa,  Canada,  18-20  October  1999,  and  published  in  RTO  MP-36. 
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1 INTRODUCTION 

Par  definition  l'analyse  aeroelastique  couple  2 types 
de  modele  qui  au  depart  sont  issus  de  calculs  sur 
plan : 

- modele  structural  dynamique  Elements  Finis  de 
1'avion  (deformations  "souples"  de  l'avion  en 
fonction  charges  ; voir  planche  1), 

- modele  aerodynamique  theorique  "linearise" 
par  methode  de  singularite  et/ou  Elements  Finis 
Euler, stationnaire  et  instationnaire  (charges 
aerodynamiques  en  fonctions  des  mouvements 
"rigides"  et  "souples"), 

qu'elle  reunit  (voir  planche  2)  en  des  modeles  de 
mecanique  du  vol  "avion  souple"  et  de  vibrations 
aeroelastiques,  eux-memes  support  des  calculs  de 
stabilite  (divergence,  Flutter)  et  de  reponse 
structurale  en  manoeuvre,  turbulence,  etc  ... 

Une  des  complexites  de  l'analyse  aeroelastique  est 
de  rendre  ces  modeles  theoriques  coherents  avec 
les  resultats  experimentaux  issus  : 

- des  essais  statiques  et  de  vibration  au  sol, 

- des  essais  en  soufflerie, 

- des  essais  en  vol,  de  reponse  en  manoeuvre  et 
d'excitations  vibratoires. 

En  reponse  a ce  besoin,  nous  avons  developpe 
progressivement,  depuis  la  fin  des  annees  70,  une 
panoplie  d'outils  de  recalage  du  modele 
structural  dynamique  et  de  la  partie 
stationnaire  des  charges  aerodynamiques  sur  les 
resultats  experimentaux. 

Ces  outils  sont  fondes  sur  une  technique 
mathematique  originale  d'identification  de  modele 
(voir  reference  2 et  3) ; ils  sont  integres  dans  notre 
logiciel  general  d'analyse  aeroelastique  ELFINI 
(voir  reference  1).  Cette  approche  est  exposee  dans 
les  § 2 et  3. 

Nous  pressentions  que  la  meme  technique  pouvait 
s'appliquer  & I'identification  des  forces 
aerodynamiques  instationnaires  et  permettre  le 
recalage  du  calcul  des  vitesses  critiques  de  Flutter 
par  des  essais  en  vol  & des  vitesses  bien  inferieures 
(voir  reference  3). 

La  demonstration  de  validite  de  cette  approche  par 
essais  en  vol  etant  evidemment  trop  risquee,  nous 
montrons  ici  (§  4)  comment  nous  avons  profite 
pour  la  tester  de  la  campagne  de  validation  des 
modeles  de  flutter  en  presence  de  jeux  mecaniques 
par  des  essais  d'une  maquette  aeroelastique 
d'empennage  "type  Airbus"  dans  la  soufflerie  S2  de 


l'ONERA  (voir  notre  presentation  sur  ce  sujet  dans 
ce  meeting,  reference  4). 

2 PRINCIPES  GENERAUX  DES  TECH- 
NIQUES MATHEMATIQUES 

D'IDENTIFI-CATION  DE  MODELE 

Nous  nous  plagons  dans  le  cas  general  ou  nous 
disposons  des  resultats  de  mesures  experimentales 
rassembles  dans  un  vecteur  ames,  et  des  memes 
quantites  c(A.),  calculees  par  un  modele 
numerique  dependant  de  parametres 
d'ajustement  X. 

L'approche  habituelle  en  identification  de  modele, 
dont  nous  allons  montrer  les  inconvenients  en 
particulier  pour  l'aeroelasticite,  est  la  methode  des 
moindres  carres,  avec  le  principe  de  minimisation 
d'une  fonction  d'erreur  E(A.)  de  la  forme  : 

E (A.)  = £ ni  (<7|  (A.)  - cTj  mes)2 

(Ili  = facteur  de  ponderation  des  differentes 
mesures) 

Cette  minimisation  est  traitee  par  une  methode 
type  Newton,  ou  la  matrice  Hessienne  de  E est 
approximee  par : 

[H]  = Zni  (Seri  / 5A.j)t  ((5ai  / 5A.j) 

dite  matrice  de  Fisher  (voir  reference  4). 

La  methode  ne  fait  intervenir  comme  calculs  lourds 
que  ceux  des  a(X)  et  des  gradients  da  / dX  qui 
sont  calcules  de  fa?on  standard  par  ELFINI  (voir 
reference  1) : 

Le  choix  des  facteurs  de  ponderation  ni  pouvant 
poser  probleme,  est  souvent  resolu  avec  la  variante 
dite  "maximum  de  vraisemblance"  (voir  reference 
5 et  6)  ou  on  itere  la  methode  des  moindres  carres 
en  prenant  comme  facteur  de  ponderation  a chaque 
iteration  l'inverse  du  carre  de  l'erreur  de 
reconstitution  de  l'iteration  precedente. 

La  difficulte  de  ces  techniques  de  moindres  carres 
est  leur  impossibilite  de  traiter  des  parametres  ou 
combinaison  de  parametres  d'ajustement  mal 
observables  par  les  mesures  (entrainant  la 
singularite  de  la  matrice  de  Fisher),  ce  qui  est  le 
cas  presque  general  avec  les  parametres 
d’ajustement  "naturels"  des  modeles 
aerodynamiques  (voir  § 3 et  4).  La  solution, 
souvent  preconisee,  de  reduire  1'espace  des 
parametres  d'ajustement  (sous  espace  des  vecteurs 
propres  a valeurs  propres  non  nulles  de  la  matrice 
de  Fisher),  n'est  pas  "physique"  car  elle  fait 
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dependre  l'espace  des  parametres  d'ajustement  des 
points  de  mesures  choisies. 

Cet  inconvenient  nous  a amene  a reformuler 
completement  le  probleme  d'identification  sous 
la  forme : 

- minimiser  Z = Z (Xj  - Xj  th)2 

Z represente  ime  distance  entre  les  valeurs 
recherchees  Xj  des  parametres  d'ajustement  et 
leurs  valeur  Xj  th,  theoriques  si  elles  sont  issues 
de  calculs  ou  presumes  s'il  s'agit  de  donnees 
physiques. 

- En  satisfaisant  les  inequations  : 

Oj  mes  - 8i  £ CTi  (X)  ^ <J|  mes  + £i 

restitution  des  mesures  par  le  modele  a une 
precision  donnee  ei. 

Cette  minimisation  sous-contraintes  de  Z est 
traitee  par  une  technique  d'optimisation 
quadratique,  ou  comme  pour  les  moindres 
carres  les  seuls  calculs  "lourds"  sont  ceux  de 
ct(X)  et  do  / 3X. 

Les  principaux  avantages  compares  aux  moindres 
carres  sont : 

- pas  de  difficult^  avec  les  parametres  mal 
observables,  ils  sont  automatiquement  kisses  a 
leur  valeur  theorique  ; on  n'est  done  pas  limite 
par  le  nombre  des  parametres  d'ajustements, 

- les  facteurs  de  ponderation  des  mesures,  ni 
plus  ou  moins  subjectifs,  sont  remplaces  par 

une  precision  de  reconstitution  des  mesures  Ei, 
qui  peut  etre  estimee  objectivement  a partir  de 
la  precision  des  mesures  et  de  celle  du  modele, 

- si  le  biais  du  modele  est  trop  grand, 
l'algorithme  donne  un  diagnostique  clair  des 
mesures  que  le  modele  n'arrive  pas  a 
reconstituer. 


3 APPLICATION  A ^IDENTIFICATION 
DES  CHARGES  AERODYNAMIQUES 
STATIONNAIRES  A PARTIR  DES 
MESURES  EN  VOL 

La  methode  a ete  appliquee  des  la  fin  des  annees 
70  pour  le  MIRAGE  2000  (voir  reference  2)  elle  a 
ete  perfectionnee  et  utilisee  systematiquement  pour 
le  programme  RAFALE  : 

- Les  parametres  d'ajustement  sont  la 
directement  les  composantes  des  champs  de 


pression  discretisees  dans  la  base  de  charges 
ELFINI  (plusieurs  centaines  de  composantes 
pour  un  champ),  pour  les  effets  de  chaque  degre 
de  liberte  de  mouvement  rigide  (incidence, 
derapage, ...,  braquages  gouvemes, ...). 

- Les  mesures  sont  directement  les  reponses  de 
jauges  de  contraintes  reparties  sur  l'avion 
(quelques  centaines),  mesurees  pendant  un 
ensemble  de  manoeuvres  d'oscillation  de 
tangage,  roulis,  derapage  a frequences 
variables,  dites  "stimuli"  (exemple  planche 
2),  qui  rendent  observables  les  effets  de  chaque 
degry  de  liberte  du  mouvement.  Pour  observer 
les  effets  non  linyaires  de  grandes  incidences  on 
complete  par  des  montees  lentes  en  facteur  de 
charge. 

Le  modele  aeroelastique  ELFINI  (voir 
planche  1)  foumit  directement  la  reconstitution 
theorique  de  ces  mesures  cr(X)  en  fonction  des 
parametres  d'ajustement  X des  champs  de 
pression  (qui  coincident  directement  avec  les 
composantes  du  champ  de  pression  discretise), 
ainsi  que  leurs  gradients  do  / 3X. 

Ces  operateurs  sont  non  lineaires  par  l'effet 
"aeroelastique"  des  modes  souples  (ces  effets 
aerodynamiques  des  modes  souples  ne  sont  pas 
ajustes  du  fait  de  leur  faible  observability  dans 
les  manoeuvres  considerees). 

Nous  montrons  des  exemples  extraits  des 
travaux  de  calibration  des  charges  du  RAFALE. 

- Un  exemple  de  manoeuvre  de  stimulus  par 
excitation  du  canard  avec,  planche  3,  la  ryponse 
ou  incidence  de  l'avion,  la  reponse  de  la  jauge 
d'emplanture  du  canard  avant  et  apres  recalage, 
le  recalage  est  la  presque  parfait  puisque  cet 
essai  a servi  de  base  au  recalage. 

- La  comparaison  essais-calculs  avant  et  apres 
recalage,  avec  des  manoeuvres  contrees  de 
roulis  complexes,  n'ayant  pas  servi  au  recalage, 
pour  des  jauges  de  derive  et  d'elevon  exteme, 
planche  4,  le  resultat  est  tr£s  satisfaisant. 

Un  point  important  a signaler  pour  l'application  de 
la  methode  est  la  necessity  absolue  de  recalage 
prealable  du  modele  Elements  Finis  statique  de 
l'avion  (voir  planche  1)  pour  garantir  la  validite  de 
l'operateur  donnant  les  contraintes  elastiques  en 
fonction  des  charges  appliques,  au  travers  duquel 
sont  observes  les  charges  aerodynamiques.  Cette 
calibration,  ou  souvent  simple  validation,  se  fait  a 
partir  des  correlations  calcul/mesures  des  reponses 
des  jauges  de  contrainte  sous  quelques  dizaines  de 
changements  d'etalonnage  au  sol  (voir  references  2 
et  3). 
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4 IDENTIFICATION  FORCES  AERODYNA- 
MIQUES  INSTATIONNAIRES 

Pour  evaluer  la  capacity  des  techniques  precedentes 
a calibrer  les  modeles  de  calcul  des  vitesses 
critiques  de  Flutter,  nous  avons  profite  de  la 
disponibilite  du  jeu  complet  de  calculs  et  de 
mesures  experimentales  realist  avec  la  campagne 
d'essais  dans  la  soufflerie  S2  de  l'ONERA  d'une 
maquette  dynamique  d'empennage  type  Airbus, 
destine  initialement  a la  validation  des  calculs  de 
Flutter  en  presence  de  jeux  mecaniques  (voir 
planche  5). 

- Parametres  d'ajustement 

En  considerant  l'equation  du  Flutter  exprimee 
dans  la  base  des  modes  propres  : 

[ (K)  - o2  (M)  - Vi  pV2  [A(cd/V)]  x = f 

et  sachant  que  nous  ne  sommes  pas  contraints 
par  leur  observability,  nous  avons  pris 
directement  comme  parametres  d'ajustement 
tous  les  termes  des  matrices  complexes 
[A(o/V)]  sous  la  forme  rationalisee  de  Karpel 
(voir  references  6 et  7). 

Les  matrices  (K)  et  (M)  sont  considerees  Id 
comme  constantes  (modes  propres  recales  avec 
les  essais  de  vibration  sans  vent). 

- Mesures 

Le  recalage  se  fait  directement  a partir  de  la 
"mesure"  en  soufflerie  des  frequences  et 
amortissements  des  seuls  2 modes  intervenant 
dans  le  Flutter  (flexion,  empennage,  rotation 
gouveme),  a quelques  points  de  pressions 
dynamiques  (1  et  2 dans  l'exemple  que  nous 
presentons). 

A noter  que  la  "mesure"  de  ces 
frequences/amortissement  resulte  de 
l'identification  par  la  methode  de  Prony  des 
poles  des  fonctions  de  transferts  entre  l'effort  de 
la  palette  d'excitation  et  les  mesures 
accelerometriques  (voir  conclusion,  develop- 
pement). 

Nous  presentons  des  resultats  caracteristiques 
menes  a partir  du  modele  aerodynamique  par 
methode  de  doublets  (maillage  planche  5) 
n'ayant  pas  subit  de  recalage  en  particulier  sans 
le  recalage  par  les  mesures  de  pressions 
instationnaires  presente  dans  la  reference  7). 

Sont  presentes  planche  6 : 


- Comparaison  calcul-cssais  sans  recalage, 
pression  critique  calculee  128000  Pa  contre 
85000  Pa  mesuree. 

- Comparaison  avec  1 seul  point  en  Pi  de 
recalage  a 38700  Pa,  la  pression  critique 
calculee  devient  73000  Pa. 

- Comparaison  avec  2 points  en  Pi  de  recalage 
(37800  Pa  et  60000  Pa),  la  pression  critique 
calculee  est  de  84500  Pa. 

A noter  cependant  que  pour  que  l'algorithme 
trouve  ces  solutions  nous  avons  yty  conduits  a 
relacher  la  precision  de  reconstitution  de 
l'amortissement  du  mode  tres  amorti  de  torsion 
voir  axe  de  developpement  (§5). 

D'autres  tests  menes  a partir  de  calculs 
aerodynamiques  volontairement  errones  (defaut 
sur  la  forme  geometrique)  donnent  des  resultats 
semblables,  ce  qui  montrerait  une  certaine 
robustesse  de  la  methode. 


5 CONCLUSION,  AXES  DE  DEVELOPPE- 
MENT 

On  imagine  sans  peine  l'interet  potentiel  de  la 
methode  de  recalage  des  forces  aerodynamiques 
instationnaires  presentees,  en  particulier  pour 
assurer  la  security  de  1'ouverture  des  domaines  de 
vol,  cela  d'autant  plus  que  la  mise  en  oeuvre  est 
relativement  simple  et  peut  couteuse  (le  calcul  des 
gradients  des  frequences/amortissement  est 
insignifiant  quand  on  travaille  en  base  reduite,  voir 
reference  1). 

On  peut  cependant  se  poser  la  question  de  la 
validity  et  la  coherence  physique  et  mathematique 
de  tels  recalages  "brutaux"  sur  les  forces 
aerodynamiques  instationnaires  generalisees.  Ce 
qui  nous  mene  a etudier  des  sophistications 
notables  de  la  methode,  principalement : 

- Prendre  les  parametres  d'ajustement  du  modele 
aerodynamique  instationnaire  directement  "a  la 
source",  en  se  rapprochant  des  coefficients  de  la 
matrice  d'influence  aerodynamique 

[SCp  / da  (<b/V)] 

reliant  les  coefficients  de  pression  Cp  sur  la 
surface  avion  aux  incidences  locales  a de  cette 
surface. 

Un  des  avantages  est  de  forcer  la  coherence 
entre  les  effets  aerodynamiques  non 
independants  (exemple  : incidence  et  vitesse  de 
pompage  avion). 
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Prendre  aussi  les  mesures  plus  "a  la  source”, 
avec  les  fonctions  de  transfert  mouvements- 
gouvemes-reponses  accelerometriques  voir 
jauges  de  contrainte,  ce  qui  a l'avantage  de 
s'affranchir  de  l'imprecision  6ventuelle  de 
l'identification  des  frequences  amortissement 
par  la  methode  de  Prony,  et  surtout  des 
difficultes  possibles  d'appariage  modes 
calcules  - modes  mesures  dans  les  situations 
complexes.  On  pourrait,  & la  limite,  exploiter 
directement  les  mesures  en  temps. 


- Integrer  en  une  seule  approche  l'identification 
des  charges  stationnaires  presentees  §3  et  celle 
des  charges  instationnaires. 

- A la  limite  integrer  aussi  la  calibration  du 
modele  structural  dynamique,  a partir  de 
stimuli  a basse  vitesse.  L'interet  serait  de 
simplifier  la  procedure,  le  recalage  du  modele 
dynamique  se  faisant  souvent*  par  des  essais 
d'excitation  par  les  gouvemes  l'avion  reposant 
sur  son  train,  ce  qui  peut  necessiter  des 
corrections  de  modele. 

* Sur  les  MIRAGE  2000  et  RAFALE  qui 
disposent  d'une  puissance  et  d'une  bande 
passante  de  servo-commande  importante. 
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PLANCHE 1 

MODELE  ELEMENTS  FINIS  POUR  CALCULS  AEROELASTIQUES 


DEFORMEES  AEROELASTIQUES 


PLANCHE  2 


ELFINI  - AEROELASTICITY 
ENSEMBLE  DES  MODELES 


* Cette  manoeuvre  nfa  pas  6t6  utilisie  pour  le  recalage 


PLANCHE  5 


MAQUETTE  AEROELASTIQUE  D'EMPENNAGE  TYPE  "AIRBUS" 
ESSAYEE  DANS  LA  SOUFFLERIE  S2  DE  L’ONERA 


Module  Structural 
Elements  Rnis 


Module  ASrodynamique 
(doublets) 


palette  d’excitation 


encastrement 


dispositif  de  securite 


gouverne 

axe 

gouverne 


PLANCHE  6 

CONFRONTATION  CALCUL  - ESSAIS  - RECALAGE 
(Frequences  - Amortissement) 


Comparaison  essais  calculs  Recalage  sur  un  point  Recalage  sur  deux  points 

avant  recalage  (flexion)  de  mesure  de  mesure 


